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Abstract 

We propose an approach to the nonvanishing boundary value problem for integrable 
hierarchies based on the dressing method. Then we apply the method to the AKNS 
hierarchy. The solutions are found by introducing appropriate vertex operators that 
takes into account the boundary conditions. 



1 Introduction 



The structure of affine Kac- Moody algebras has provided considerable insight underlining 
the construction of integrable hierarchies in 1 + 1 spacetime dimensions. An integrable 
hierarchy consists of a series of equations of motion with a common algebraic structure. 
They are constructed out of a zero curvature condition which involves a pair of gauge 
potentials, lying within the affine Lie algebra. A crucial ingredient is the decomposition 
of the affine Lie algebra Q into graded subspaccs, determined by a grading operator 
Q. The integrable hierarchy is further specified by a constant, grade one semi-simple 
element E^^\ The choice of Q, Q and E^^^ are the basic ingredients that fully determine 
and classify the entire hierarchy (see for instance ref. [1] for a review). 

Representation theory of Kac-Moody algebras play also an important role in con- 
structing systematically soliton solutions of integrable hierarchies. The main idea of the 
dressing method consists in, by gauge transformation, to map a simple vacuum con- 
figuration into a non trivial one- or multi-soliton solution. The pure gauge solution of 
the zero curvature representation leads to explicit spacetime dependence for the vacuum 
configuration which in turn, generates by gauge transformation, the non-trivial soliton 
solutions of the hierarchy [2, 3]. It becomes clear that the sohitions are classified into 
conjugacy classes according to the choice of vacuum configuration. In general, the vac- 
uum is taken as the zero field configuration and the soliton solutions are constructed and 
classified in terms of vertex operators [4] . 

A dressing method approach to nontrivial constant vacuum configuration was pro- 
posed in [5] when considering the negative even fiows of the modified Korteweg-de Vries 
(mKdV) hierarchy, which do not admit trivial vacuum configuration, i.e. vanishing 
boundary condition. A deformed vertex operator was introduced and the method was 
further applied to the whole mKdV hierarchy with nonvanishing boundary condition 
and to a hierarchy containing the Gardner equation [6]. Here in these notes, we shall 
provide further clarifications discussing nonvanishing boundary conditions for integrable 
hierarchies in general, constructed with the structure proposed in [7] (and references 
therein). We point out that not all the models within the hierarchy admit solutions 
with nonvanishing (constant) boundary condition. We then consider as an example the 
Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy with nontrivial constant vacuum config- 
uration. Results for the focusing nonlinear Schroedinger (NLS) equation were recently 
obtained in ref. [8]. A hybrid of the dressing and Hirota method for a multi-component 
generalization of the AKNS hierarchy was also considered recently in [9]. 

We first describe the general aspects of the dressing method in section 3, together 
with a discussion on the possible different boundary conditions. We point out the re- 
quirements for the individual models of the hierarchy to admit a nonvanishing constant 
boundary value solutions. Later in section 4, we discuss the construction of the AKNS 
hierarchy, followed by the construction of its solutions, section 5. In section 6 we present 
the deformed vertex operators used to construct solitonic solutions with nonvanishing 
boundary condition and illustrate with explicit examples. 
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2 Algebraic concepts and integrable hierarchies 



In this section wc introduce the main algebraic concepts used for the construction of 
integrable models and to obtain its solutions. A more detailed exposition can be found 
in [10, 11]. 

Let ^ be a semi-simple finite-dimensional Lie algebra. The infinite-dimensional loop 
algebra L{Q) is defined as the tensor product of Q with integer powers of the so called 
complex spectral parameter A, L{g) = g^C[X, X'^]. If e ^ then = Ta<^X^ G L{g) 
for n e Z. The commutator of the loop algebra is given by 

lT-,T,^] = [r„,r,]®A"+-, (1) 

where [T(i,T;,] is the Q commutator. 

The central extension is performed by the introduction of the operator c, which 
commutes with all the others [c, T"] = 0. Furthermore, consider the spectral derivative 
operator d = such that [c?, T^] = n-T^, so it measures the power of the spectral 
parameter. The affine Kac-Moody Lie algebra is defined by ^ = L{g) ® Cc© Cc? and for 
T", e g the commutator is provided by 

[t:, T^] = [Ta, n] ® A'^+'" + n6n+„,fl{Ta\n)c, (2) 

where {Ta\Tb) is the Killing form on g. 

The algebra g can be decomposed into Z-graded subspaces by the introduction of a 
grading operator Q such that 

G = ^G^'\ (3) 

where g^'^ = {T^ G g \ [Q,T^] = jT^}. The integer j is the grade of the operators 
defined with respect to Q. It follows from the Jacobi identity that if T" e g^^^ and 
e 0^^^ then [7^,7^] G g^'+^\ 

Let be a semi-s imple element of g. Its kernel is defined by /C = {T^ eg\[E, T^] = 
O}. Its complement is called the image subspace A4 and g = JC (B M- It can be readily 
verified from this definition and the Jacobi identity that [/C,/C] C /C and [/C, C A4. 
We also assume the symmetric space structure [A^, A^] C /C. 

Consider the linear system in 1+1 spacetime dimension 

(9, + C/)* = 0, {dt + Vp^O, (4) 

where U,V G g and ^ is an element of the Lie group of g. The compatibility of this 
system ensures integrability, yielding the zero curvature equation 

[d, + U,dt + V]=0. (5) 

The general structure to construct nontrivial integrable models through an algebraic 
approach was discussed in [7, 10] (see also references therein). Let Q be a grading 
operator and choose a constant, independent of x and t, grade one semi-simple element 
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E = E^^\ Let A^^\x, t) be a linear combination of zero grade operators lying entirely in 
M.. The coefficients of this linear combination are spacetime functions, i.e. the fields of 
the model. Therefore, if U and V have the following algebraic structure [7] 

U = E^^^ + A^''\x,t), (6a) 
V=J2 i^^t)^ D^"^^ + + • • • + D^-P^ , (6b) 

i=—p 

where p,q > arc nonnegative integers, the zero curvature equation (5) can be solved 
grade by grade, determining each D^^^ in terms of the fields in A^^\x,t); the zero grade 
projection will provide the nonlinear integrable field equations. The explicit examples 
considered in the following sections will illustrate this assertion. The particular case 
where V = D^~^^ is equivalent to a relativistic model obtained from a Hamiltonian 
reduction of the Wess-Zumino-Novikov-Witten model. The infinite number differential 
equations obtained from this algebraic structure, one for each choice of p and q, constitute 
an integrable hierarchy immediately classified hy Q, Q and E^^\ 

Consider a unitary highest weight representation of G, where the hermitian relations 
are (T")^ = T"". There exists a finite set of states {|/ifc)}, A; = 0, 1, ... , rank^, that are 
annihilated by every operator with a positive power on the spectral parameter, 

T:\f,,) = 0, n > 0. (7) 

The fundamental state I/Iq) is annihilated by every operator, including {T'a }, except the 
central term. Under the action of c all states obey 

c\fj.k) = l/Ufe)- (8) 



3 Dressing method 

Let us briefly review the dressing method [2, 3]. Assume a known solution of the zero 
curvature equation (5), in pure gauge form, called vacuum 

C/o = -9,*o*o\ l^o = -5**0*0 (9) 
The dressing method arises as a factorization problem 

QZ'e+ = *o^?*o \ (10) 

where g is an arbitrary constant group element, independent of x and t, and 0± are 
factorized into positive and negative grade operators, respectively, according to a grading 
operator (3), 

©_ — exp(m*-~^-') exp(m''~^'') • • • , B+ = exp(m''°-') exp(m''^-') • • • . (11) 

where m^^^ e Q^^K The factorization problem (10) ensures that, the dressing operators 
Q± map a vacuum solution into a new solution by the gauge transformations 

u = e^UoQi' - d^e±e^\ (m) 

V = e±VoQl^ - dtQ±ei\ (12b) 
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3.1 Vanishing boundary condition 

Taking a zero field configuration A^^^ ^ in (6) we have the constant vacuum operators 



Uo^E^^^ = n^^\ (13a) 
Vo = dS^^ + Dj-"^ = + fl^-P') . (13b) 

where il^^^ G /C. This conclusion comes from the projection of the zero curvature equation 
into each graded subspace and the resulting equations for the highest and lowest grades 
[7]. Using (13a) in the gauge transformation (12a) with 0_|_, taking into account the 
central extension^ of the loop algebra denoted by — i^^c, and projecting the resulting 
equation into a zero graded subspace we get 

- u^c = -a^exp(m^°)) exp(-m^°)). (14) 

Assume a pure gauge connection, involving zero grade generators only, for the field 
operator 

= -B^B-\ (15) 

Substituting (15) into (14) gives 

exp (m'-"-') = i? cxp (z/c) . (16) 

Consider the fundamental state |yUo) and generic states luk), l/J-i) that are not annihilated 
by some operators T^. In fact, the generic states are chosen as the ones that survive 
under the action of the zero grade operators associated with the fields in A^^\ This 
will be clear when we consider a concrete example. Therefore, projecting (10) between 
appropriate highest weight states, we get the implicit solution of the model 

{tik\B\i,i) = ^, (17) 
Too 

where we have defined the tau functions 

Too = (/iol^O^^O Vo), Tfc/ = (/ifcl^o^^'o Vi)- (18) 

The vacuum operators (13) will commute, since for zero field configuration they lie 
in /C, so integrating (9) yields 

*o = exp{-Q(^)x - + n^-P^)t}. (19) 

Soliton solutions arise by choosing g in terms of vertex operators in the form [4] 

N 



g = l[exp{X,). (20) 



1=1 



^In the construction of integrable models only the loop algebra plays a significant role, but to employ 
the highest weight representation the central term needs to be carefully considered. The equations of 
motion are invariant under U U + fc. We are choosing / = —v^ for a convenient notation. 
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Suppose the vertex is an eigenstate of the vacuum, i.e. 

= -r]i{x,t)Xi. (21) 

The function rji contains the dispersion relation in an exphcit spacetime dependence. 
Then, the tau functions (18) can be computed in an exphcit way, 

(N \ N N 

n '''' *o ' = n ^oe^^^o ' = n (22) 
i=l ) i=\ i=\ 

Using the expansion e^Be"^ — S + [A, S] + [A, 5]] + ••• with the eigenvalue 
equation (21) we obtain the A^-soliton solution 

AT N 

i=i i=i (^23) 

+ ai2e''i+^2 + ai3e^i+''3 + a23e^2+''3 + 01236^1 +"^+^3 + • • • , 

where we have assumed the nilpotency property of the vertex, (Xj) "l/x;) = for 
n >2, and defined the following matrix elements 

Gi = {nklXilni), Uij ^ {iJ,k\XiXj\fj,i), Gijk = {iJ,k\XiXjXk\ni), etc. (24) 
3.2 Nonvanishing boundary condition 

Let us now assume a constant, nonzero, field configuration A^^^ — > Aq^^ Then (6) yields 
the constant vacuum operators 

[/o = ^«+4°\ (25a) 

Vo = D^^^ + Di''-'^ + --- . (25b) 

The analogous of eq. (14) becomes 

- u^c = exp(m(°))4°^ exp(-m(°)) - exp(m(°)) exp(-m(°)) . (26) 

Assuming the form (16), then (26) yields a relation between the physical fields in A^^\ 
functionals associated with zero grade generators in B and the nonzero vacuum, i.e. 

= BA^^^B-^ - B^B-\ (27) 

When ^ we recover the relation (15). The general solution of the model is still 
given in the form 









(/xol^'o^^'o'l 


^0) 
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but now with ^'q satisfying the hnear system (9) with the more general vacuum operators 
(25) and henceforth, depending upon the constant vacuum solution. We emphasize that 
the implicit solution (28) and (27) is valid for an arbitrary vacuum solution, even if it 
is a general spacetime function. For the general case, however, the difficulty arises in 
integrating the linear system (9), while for a constant vacuum solution this is trivial. 

Assuming a constant vacuum, the operators (25) does not lie purely in /C anymore, 
but they are supposed to be solution of the zero curvature equation, i.e. [Uq, Vq] = 0, so 
integrating (9) gives 

^o^exp{-Uox-Vot). (29) 

Assuming again that g has the form (20) and that there exists vertex operators satisfying 
eigenvalue equations like (21), but now with operators (25), it will be possible to compute 
the tau functions in (28). Because (25) depend on the constant vacuum fields it suggests 
that for (21) to be valid, Xi will also depend on these constants fields, so the dispersion 
relation and multi-solitonic interaction terms will depend on the boundary conditions. 

3.3 Models admitting nonvanishing boundary value solutions 

Not all the models in the hierarchy defined by (6) have solutions with nonvanishing 
boundary condition. The vacuum must be a solution of the zero curvature equation, 
but for this to be true, in the case of constant boundary condition, a restriction on the 
structure of the Lax pair is required. 

The operator (6b) depends on the fields. When taking the fields to vanish, only 
the operator with highest (and/or lowest) grade that is a constant operator purely in 
JC remains, see (13b). Clearly, in this situation there is no restriction besides the own 
algebraic construction (6). 

For a nonzero constant field configuration, other operators in (25b) that previously 
vanished now remain, and the whole combination should commute with (25a). In (25a) 
we have the structure il^^^ = )C^^^ + If (25b) keeps the same algebraic structure, 

i.e. if it is formed by a linear combination of operators of the kind 

qU) = jci^) + >jO-i) = X^i-^ni^) (30) 

where j < n and A is the loop algebra spectral parameter, then we will have [Uq, Vq] = 0. 
Note that in (30) it is required at least two operators differing by one grade. This will not 
be possible for the first model of the negative hierarchy, V = 'Yldl-p ^^at contains 
only one operator V = D^~^\ These are the relativistic models of the hierarchy, like 
sinh-Gordon and Lund-Regge for instance. In sum, at least the relativistic models of any 
hierarchy constructed with the structure (6) will not have solutions with nonvanishing 
(constant) boundary condition. 

4 The AKNS hierarchy 

The previously discussed approach for a nontrivial vacuum solution was considered for 
the mKdV hierarchy in [5, 6]. Now we present another example by considering the non- 
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abelian AKNS hierarchy. Let Q = Ai ^ s^2 with the following commutation relations 

[H^, H^] = 2n5„+„,o£, [i^", E^^] = ±2^2+-, [K^ ^-a] = H^'^"' + n5n+m,oC (31) 

Consider the homogeneous grading operator Q = d, which decomposes the algebra into 
the graded subspaces 

gU) ^ ^H^^ El Ei^]. (32) 

Choosing the constant semi-simple element in (6a) as E^^'> — H^, we have A^^'> = qE^ + 
rE^^ where q — q{x, t) and r = r{x, t) are both fields of the models within the hierarchy. 
Let 

D^^'^ ^ajH^ + bjEi + cjEi^ G Q^^^ (33) 

where functions of x and t. Neglecting for the moment the central terms (i.e. 

c = 0) and considering only positive grade operators in (6b) we have the zero curvature 
equation for the positive AKNS hierarchy 

[d, + H^ + qEl + r£;°,, 9* + D^") + D^""^) + • • • + L>(°)] = 0, (34) 

where n = 1, 2, 3, ... is an arbitrary positive integer. Solving (34) grade by grade, starting 
from the highest one, we can determine all coefficients in (33) in terms of the fields q 
and r. The zero grade projection provides the field equations. 

For the case where n = 2, fixing the two arbitrary constant coefficients as 02 = 1 and 
ai = 0, and ignoring possible integration constants as usual, we get the famous AKNS 
system 

Qt = -Iqxx + q'^r, n = Ir^x - qr'^ ■ (35) 
Choosing q = ip, r = cr^p*, t — )■ —it and x ^ ix we obtain the well known NLS equation 

m-liPxx-cr\(pf^p^O. (36) 

The cases a — ±1 correspond to the focusing and defocusing NLS equations, respectively. 
Note that the nonzero constant fields q ^ qo, r ^ ro are not solution of (35). If we do 
not ignore an integration constant arising when solving for the ao coefficient in (33) , but 
denote it conveniently as ^qoVo, we get the system 

qt = -^qxx + {qr - qoro)q, n = |r^x - {qr - qoro)r, (37) 
whose Lax pair is given by 

U^H^ + qE^^ + rE''_^, (38a) 
V^H' + qEl + tEI^ -\{qT- qoro)H' - ^q^E^ + ir,£;°,. (38b) 

Note that now q ^ qo, r ^ tq is a solution of (37). The choice q = Lp, r = aLp*, t — )■ —it 
and X ^ ix yields the focusing and defocusing NLS equations admitting nonvanishing 
(constant) boundary value solution 

m - \^xx - (T(|<^r - A^) = 0, (39) 
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where = |(/?oP is a real constant. Solutions of this equation with a = 1 (focusing) was 
the object under study in the recent paper of Zakharov and Gelash [8] . 
Solving (34) for the n — 3 case we get 

qt = Iqxxx - Iqrqx, n = ^r^xx - Iqrr^, (40) 

whose Lax pair given by (38a) and 

V = H^ + qEl + rE\ - \qrH^ - \q^El + \r^E\ 

+ - ?^-)^° + - 2?'0^a + - 2q\)E^_^. (41) 

The constant fields q ^ q^^ r ^ are solution of (40). Note that taking r = ±g in 
(40) yields the niKdV equation, r = ±1 the KdV equation and r — a -\- fiq the Gardner 
equation with arbitrary coefficients 

4% = qxxx - ^aqq^ - Qpq^qx- (42) 

Solutions with nonvanishing boundary condition for the mKdV hierarchy and solutions 
of a deformed hierarchy containing the Gardner equation were considered in [6] . 

The first negative flow of the AKNS hierarchy is the Lund-Regge model, that as 
discussed in section 3.3, does not have nonvanishing (constant) boundary value solutions. 



5 AKNS solutions 

Let 

=H^^ ^^^n-l ^ (43) 

Taking q ^ q^, r ^ Tq in (38) and (41) yield the respective vacuum Lax operators 

Uo = n^^^ for (38a), (44a) 

Vo = for (38b), (44b) 

Vo = n^^^ - yoro^^^^ for (41). (44c) 

Let 

B = exp(xi?° J exp(0i/°) exp(^£;°). (45) 
Solving (27) will imply in the following relations 

0. = ro^(l + x^e""^) - qoxe^'^ + X^xC^t 
g = e2*(^o-ro^'-^x), (46) 
r = roe-^t' + ^'''{q^x^ - rox'^') - Xx - xV.e'^ 

Introducing the auxiliary functions 

X = Xe^ ^ = M A=l + x^, (47) 
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the expressions (46) can be further simphfied to 

g=|e^'^-fe^ r = roAe-2<^ - x.e-^. (48) 

Taking into account the central terms (i.e. c 7^ 0), consider the highest weight states of 
5^2; {|a*o)) |a*i)}- These states obey the following actions, 

H^W) = 0, ElJ^ik) = 0, (forn>0) (49) 

E^Jfik)^0, i^Vo)=0, i/Vi) = lA^i), c\iik) = W)- (50) 
The adjoint relations are 

{H^y = H-^, {El,y = E-2, c^^c, S = d. (51) 
Besides the highest weight states we use the same notation to denote the particular state 

= E'_M + 0. (52) 
Using these states in (28) with (45), we have the following tau functions 



Too = e"" = (^ol^ofi-^o Vo), 
rii = e'^+^ = (/xi|^o<7^oVi), 

T2i = xe'^+* = (/X2|*o^*o"Vi), 



(53) 



from which we have 



e* = I^. x = ^. ^ = ^, X = ^, ^ = ^. (54) 

Too Til Til Too Too 

Substituting the tau functions in (48) we get the explicit dependence for the nonvanishing 
boundary value solutions of the AKNS hierarchy 

, = *^^ + ^<5^1-?p^, (55a) 

'''00 + T12T21 Too Tqq + Ti2r2i 

r = ro"°°+:""" + (55b) 

Til TooTii 



6 Vertex operators 

We now introduce vertex operators that are eigcnstates of the nontrivial vacuum Lax 
operators (44). Consider g in the usual solitonic form (20), with vertex operators defined 
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by 



where X- 



(n) 



X, 



(n) 



gon) 



(n) 



where K' 



(n) 



K. 



n ro{Kj-qoro) 



Ki[Kf + qoro) Ki 



(56a) 
(56b) 
(56c) 
(56d) 



where Ki is a complex parameter. In the hmit go, ro — >■ these vertex operators recover 
the well known vertices used for zero boundary condition [12]. A direct calculation will 
give the following eigenvalue equations, for Q^") defined by (43), 



p^\Xi]=- 



^1 Ki 



-\ Ki- 



gprp 

Ki 

gpyp 

Ki 



n-1 



n— 1 



Ki + 



Ki + 



gpro 


Xi, 


Ki 




gprp" 


Yi. 


Ki 





(57a) 
(57b) 



Therefore, the spacetime dependence of the solutions of system (37) will be given by 

go?^p\ 1 

1 ' 



Vi^+i Ki + -^j 



x+ - \ Ki - 



K7 



= - K , 



Ki 



X 



2_ (gp^p): 

Ki ^2 ) '^^ 



(58) 



while for the system (40) we have 



gpj'o 



Ki I 4 

> _ , , gp'^p \ 1 

?i — \ Ki-\ I X — 

Ki J A 



^ + 7 '^i ~ ^KiqoTo - 3 



(go?'o)^ , (go?^p) 



+ 



Ki 



K? 



Kf - SKiqoro - 3 



(go^^p)^ , (gp^'o) 



Ki 



+ 



Ki 



(59) 



t. 



The matrix elements involving one vertex operator only are 



{llQ\Yi\llQ) 

{lii\Yi\iii) 



ro{Ki - gprp) 
Ki (kI + gpro) ' 

-2Kiro 



Kt + gpro ■ 

(/^l|>"i|^2) = 1, 

{^^2\Yi\|ll) = — I 



(/iol^il/ip) 
{lii\Xi\iii) 

{lll\Xi\ll2) 



QoiKj - gpro) 



Ki (k • + gpro) ' 

-2goro 



Ki + gpro) 



.g| 

Kf 



K. 



2 ' 



{lJL2\Xi\lJix) = 1. 



(60a) 
(60b) 
(60c) 
(60d) 
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Unlike the case of vanishing boundary condition, the solution involving only one vertex 
is not trivial. Choosing g = e^^ yields 

Too = 1 + {fJ'o\Xi\fio)e'^\ (61a) 

rii = l + (/ii|Xi|/xi)e''\ (61b) 

ri2 = (/Ui|Xi|/i2)e''\ (61c) 

r2i = (//2|Xi|//i)e''^ (61d) 

Analogous expressions are also valid for Yi, replacing r/i — )■ ^i. Substituting the tau 
functions (61) in (55) we get the explicit 1-vertex solution of the positive AKNS hierarchy. 

For a solution involving two vertices like g = e^^e^^, we have more complicated 
matrix elements, but its general form is given by 

Too = 1 + (^ol^il^o)e''^ + (/io|r2|/io)e«^ + (/io|XlF2|/io)e''^+«^ (62a) 

rn = 1 + {fii\X^\^i^)e^^ + {i2^\Y,\i2^)e^' + {fi,\X,Y,\i2,)e'^'+^\ (62b) 

Ti2 = (//i|Xi|/i2)e''^ + {fii\Y2\fi2)e^' + {fii\XiY2\fi2)e"'^^\ (62c) 

r2i = (//2|^i|//i)e''^ + {fi2\Y2\fXi)e^' + {fi2\XiY2\fXi)e^'+^\ (62d) 

Besides the matrix elements involving one vertex operator (60), the nilpotency property 
was explicitly verified and the following two vertices matrix elements calculated, 

[Kf + go^O ) [K^j + goro j KiKj [Ki - Kj) 

{^^\X^YM) - (^°^"("--"^)'+("-"^-^°^°)'}("-"^+^°^")', (63b) 
= -2.0 , j:-"--;?!'"-"-;*'-''' (63c) 

[Kf + qoro)[K^j + qoTQ) [Ki - KjjKf 

(/X2|X,F,|/.i) = -2ro ' \ °^ (63d) 

[k; + goroj (k^ + go^oj (Kj - Kj)K^j 

The matrix elements (63) and (60) recover known results when go, '"o 0, corresponding 
to the vanishing boundary condition [12]. 



7 Conclusions 

We have shown how to employ the algebraic dressing method, using vertex operators, to 
construct nonvaninshing (constant) boundary value solutions for integrable hierarchies. 
This method was applied previously for the mKdV hierarchy [5, 6] and in this letter for 
the AKNS hierarchy, completing the possible classical models obtained using the algebra 
5^2 with the construction of integrable hierarchies [7]. We have algebraically explained 
how the dispersion relation and multi-soliton interaction depend on the boundary con- 
ditions, here arising as a consequence of the vertex operators (56). 
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Our solutions are quite general, solving at once various individual models contained 
in the AKNS hierarchy hke the focusing and defocusing NLS, mKdV, KdV and Gardner 
equations. We also have obtained an explicit 2-soliton solution, involving two vertex 

operators. Of course higher multi-soliton solutions can be constructed following the 
same approach, we only face a technical difficulty in computing the matrix elements but 
the procedure is well established. 

Acknowledgements. We thank CNPQ and Fapesp for financial suppport. We would 
like to thank L. Feher for suggestions. 
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